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ABSTRACT 

We present a new analysis of the motion of pressure-confined, broad line region (BLR) clouds in 
active galactic nuclei (ACNs) taking into account the combined influence of gravity and radiation 
pressure. We calculate cloud orbits under a large range of conditions and include the effect of column 
density variation as a function of location. The dependence of radiation pressure force on the level 
of ionization and the column density are accurately computed. The main results are: a. The mean 
cloud locations {tblr) and line widths (FWHMs) are combined in such a way that the simple virial 
mass estimate, tblrFWHM 2 /G, gives a reasonable approximation to Mbh even when radiation 
pressure force is important. The reason is that L/M rather than L is the main parameter affecting 
the planar cloud motion, b. Reproducing the mean observed Tblr> FWHM and line intensity of H/3 
and C IV AI549 requires at least two different populations of clouds, c. The cloud location is a function 
of both L 1 / 2 and L/M. Given this, we suggest a new approximation for tblr which, when inserted 
into the BH mass equation, results in a new approximation for Mbh- The new expression involves 
L 1 / 2 , FWHM and two constants that are obtained from a comparison with available M — c* mass 
estimates. It deviates only slightly from the old mass estimate at all luminosities, d. The quality 
of the present black hole mass estimators depends, critically, on the way the present M — tr* AGN 
sample (29 objects) represents the overall population, in particular the distribution of L/L^dd- 
Subject headings: Galaxies: Active - Galaxies: Black holes - Galaxies: Nuclei - Galaxies: Quasars: 
Emission Lines 



1. INTRODUCTION 

The profiles of the broad emission lines in the spectrum 
of active galactic nuclei (AGNs) are the main source of 
information about the motion of the high density gas 
in the broad line region (BLR). Detailed studies of such 
profiles have been the focus of intense investigation for 
many years (see Netzer 1990 for a review of older work 
and Marziani et al. 1996 and Richards et al. 2002 for 
more recent publications). Unfortunately, several rather 
different geometries can conspire to result in similar line 
profiles and today, there is no way to infer, directly, the 
global BLR motion from line profile fitting. 

A less ambitious goal is to use a measure of the ob- 
served line width, e.g. the line FWHM, or the line disper- 
sion (see Peterson et al. 2004 for definitions) as indica- 
tors of the mean emissivity- weighted velocity of the BLR 
gas. Such measurements are crucial for deducing black 
hole (BH) mass (Mbh) in cases where the emissivity- 
weighted radius, tblr, is measured directly from rever- 
beration mapping (RM) experiments, or estimated from 
from h-rsLR relationships that are based on such studies 
(see Kaspi et al. 2000; Kaspi et al 2005; Vestergaard and 
Peterson 2006 for reviews). A typical expression of this 
type is 

£5100 



and 7 = 0.6 ±0.1. The constant a depends on the line in 
question. For H/3, a ~ 0.4 pc (e,g, Bentz et al. 2009) and 
for CivA1549, a ~ 0.13 pc (Kaspi et al. 2007 after as- 
suming £1350 = 2^5100). For a virialized BLR, the above 
t blr can be combined with a measure of the FWHM, or 
the line dispersion, to obtain the BH mass, 



M BH = f tblrFWHM 2 /G 



(2) 



tblr 



10 46 ergs" 



pc, 



(1) 



where L5100 is the continuum luminosity (ALa) at 5100A 
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where the constant / is a geometrical correction factor of 
order unity that takes into accounts the (unknown) gas 
distribution and dynamics. Various possible values of / 
have been computed by Collin et al. (2006) for various 
possible geometries. However, the only empirical way 
to determine / is to compare the results of eqn. [5] with 
independent measurements of Mbh, like those available 
in the the case where the central BH resides in a bulge 
and Mbh can be estimates from the M — a* relationship 
(e.g. Tremaine et al. 2002). Such comparisons by Onken 
et al (2004) and by Woo et al. (2010), using the H/3 RM 
data base, suggest / = 1 ± 0.1. 

In a recent paper, Marconi et al. (2008; hereafter M08) 
investigated the role of radiation pressure force and its 
effect on the motion of the BLR gas and the required 
modification to the BH mass estimate. According to 
M08, radiation pressure plays an important role in af- 
fecting the cloud motion provided the column density 
(N co \) of most BLR clouds is smaller than about 10 23 
cm~ 2 . According to M08, in such a case, there is a need 
to add a second term to eq. [2] This term depends on the 
source luminosity and N co \. The modified form suggested 
in M08 is 



M BH = hr BL RFWHM 2 /G + f g L/N c , 



(3) 
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where/i replaces / in eqn. [5] and f g is a second constant. 
If L = L5100/IO ergs -1 and A co i is measured in units 
of 10 23 cm- 2 , f g ~ 1O 7 ' 7 M : According to M08, failing 
to account for the second term results in the underesti- 
mation of A/bh- Obviously, the inclusion of such a term 
results in /1 < /. M08 repeated the analysis of Onken et 
al. (2004) and Vestergaard and Peterson (2006), taking 
into account the new term and solving for /1 and JV co i. 
This resulted in fx ~ 0.56 and A col ~ 10 23 cm" 2 . 

The M08 suggestion can be tested by comparing low 
redshift samples of type-I and type-II AGNs since the es- 
timate of Mbh in the latter does not involve the source 
luminosity and gas dynamics. Netzer (2009) carried out 
such a comparison and found that radiation pressure 
force plays only a marginal role in such sources. The 
conclusion is that, in many AGNs, the mean column den- 
sity of the BLR clouds exceeds ~ 10 23 cm~ 2 . In a later 
work, Marconi et al (2009; hereafter M09) argued that 
firm conclusions regarding the role of radiation pressure 
force are difficult to obtain since the column density in 
some BLRs can be different than in others and there is 
no simple way to evaluate the overall effect of such a 
column density distribution. The treatment of a certain 
type of cloud in all sources, or even in a single BLR, is 
of course highly simplified and eqns. [5] and [3] must be 
treated as crude first approximations. 

The critical and detailed evaluation of the role of radi- 
ation pressure force in "real" BLRs is the subject of the 
present paper. In §2 we present our basic equations and 
in §3 we use them to calculate various expected broad 
emission line profiles and mass normalization factors, /. 
§4 deals with the evaluation of present day Mbh esti- 
mates and suggests a new way to estimate Mbh and tblr 
which is consistent with our calculations. 

2. CLOUD MOTION IN THE BLR 

In this work we focus on the "cloud model" of the 
BLR. The general framework of this model is explained 
in Netzer (1990) and in Kaspi and Netzer (1999) and 
a major empirical justification is obtained from the re- 
cent X-ray detected single blobs, or clouds, moving in 
a region which is typical, in terms of velocity and di- 
mension, of the BLR (Risaliti et al. 2010; Maiolino et 
al. (2010). We do not consider the locally optimally- 
emitting cloud (LOC) model (Baldwin et al. 1995; Ko- 
rista & Goad 2000) where, at every location, there is a 
large range in cloud properties. The dynamics of the 
BLR gas in this model has never been treated and is 
far more complicated than the one considered here. An- 
other possibility that has been discussed, extensively, is 
that wind from the inner disk plays an important role in 
feeding and driving the BLR gas. Possible evidence for 
this scenario comes from radio observations (e.g. Vester- 
gaard, Wilkes, & Barthel 2000; Jarvis & McLure 2006) 
and spectropolarimetry (Smith et al. 2004; Young et 
al. 2007). Theoretical considerations are discussed in 
Bottorff et al. (1997), Murray & Chiang (1997), Proga, 
Stone, & Kallman (2000), Everett (2003), Young et al. 
(2007), and several other papers. While our calculations 
apply to any cloud, even those created and driven by 
such winds, the specific examples given below are more 
applicable to bound clouds where inward and outward 
motions are both allowed. 



2.1. The equation of motion of BLR clouds 
The basic equation of motion, ignoring drag force, is 



a(r) 



-^[M(r)-l/r]-if , 
/im#c47rr z p dr 



(4) 



where M(r) is the force multiplier, Lm is the bolometric 
luminosity, /i is the average number of nucleons per elec- 
tron, and r =L/LEdd- The force multiplier depends on 
the gas composition and its level of ionization. An inter- 
esting case is a Compton thin neutral cloud that absorbs 
all the ionizing radiation (a Compton thin "block"). In 
this case M(r) ~ a{r)/ (<jtN co i) where N co i is the hy- 
drogen column density and a(r) is the fraction of the 
bolometric luminosity which is absorbed by the gas. For 
such a "block", a(r) = Li on /Lboi but in general a(r) is 
radius dependent because of the changing column density 
and level of ionization of the gas (see below). 
Ignoring thermal pressure we obtain 



1.14 x 10~ u a(r) 8.8x10 



-13 



N- 



2:i 



r 



(5) 



where A23 = A CO ;/10 23 . Thus, radiation pressure is the 
dominant force when 



T > 7.7 x 10~ 2 



-2^ 



23 



a(r) 



(6) 



The above expressions, including the one for the limiting 
r, include only radial terms and assume a pure radially 
dependent radiation pressure force. The calculations of 
real orbits, and the conditions for cloud escape, require 
their integration and will thus include the standard con- 
stants of motion (energy and angular momentum). Ob- 
viously, the conditions for escape depend on the cloud 
azimuthal velocity, v$, and can differ substantially from 
what is obtained by using eqn. [5] 

M08 and M09 derived similar expressions for the case of 
completely opaque clouds. According to them, radiation 
pressure dominates the cloud motion if 

r > 1.27 x 10- 2 & 5 iooA 23 (7) 

where 65100 = -k&oZ /-^5ioo ■ The two expressions provide 
the same limiting T when 

a(r)& 5 ioo -6.1 . (8) 

A recent paper by Ferland et al. (2009), where mostly 
neutral, infalling clouds are considered, reaches basically 
identical conclusions. 

2.2. Confined clouds 

BLR clouds are likely to be confined. The confin- 
ing mechanism is not known but high temperature gas 
and magnetic confinement have been proposed. The ap- 
proach chosen here is consistent with the idea of mag- 
netic confinement and some justifications for it are given 
in Rees, Netzer and Ferland (1989). We adopt a sim- 
ple model of numerous individual clouds that are mov- 
ing under the combined influence of the BH gravity and 
radiation pressure force. Following Netzer (1990), and 
Kaspi and Netzer (1999), we assume that clouds retain 
their mass as they move in or out and the gas density 
changes with radius in a way that depends on the radial 
changes of the confining pressure. 



Radiation pressure in AGN 



3 



Assume the external pressure and hence the gas den- 
sity in individual clouds are proportional to the radial 
coordinate, nu oc r~ s . A reasonable guess that agrees 
with observations is 1 < s < 5/2 (Rees, Netzer and Fer- 
land 1989). This results in a radial dependence of the 
ionization parameter (the ratio of ionizing photon den- 
sity to gas density), U oc r s ~ 2 . For spherical clouds, 
N co i oc r~ 2s / 3 , R c oc W 3 and A c oc r 2s / 3 , where R c 
is the cloud radius and A c its geometrical cross section. 
The line intensity contributed by a single cloud, e(r), de- 
pends on its covering factor and the line emissivity j(r) 
which depends on the conditions in the gas, 



e(r) oc j(r)A c /r 2 ocj(r)r 2s / 3 ~ 



(9) 



In the real calculations we ignore factors of order unity 
relating the mean cloud "size" and its mean column den- 
sity since this is not known and require different type of 
calculations. 

The above considerations suggest that the importance 
of radiation pressure increases with distance because of 
the dependence of N co i on r, i.e. 



r Hm oc r- 2s / 3 /a(r) . 



(10) 



Since T depends on the global accretion rate which has 
little to do with cloud properties, the more physical ap- 
proach is to consider the case of a certain T and follow 
the cloud motion. The examples discussed below follow 
this approach. 
In this work we consider three types of clouds: 

1. Very large column density clouds where radiation 
pressure force is negligible at all distances. Here 
virial cloud motion is a good approximation (the 
Ferland et al 2009 infalling clouds belong to this 
category) . 

2. Cloud for which radiation pressure is very impor- 
tant somewhere inside the "classical BLR" (e.g. in- 
side the RM radius). Such clouds will escape the 
system on dynamical time scales and their contri- 
bution to the line profiles is small except for times 
immediately after a large increase in L{, i- 

3. Clouds for which radiation pressure is non- 
negligible but is not strong enough to allow escape. 
Such clouds are the ones discussed by M08 albeit 
without the radial dependence of N co i considered 
here. This case is the one most relevant to real 
BLRs and we discuss it in detail in the following 
section. 

2.3. Modified equation of motion 

The modified equation of motion is obtained from 
cqn. [5] by including the radial dependence of N co i . We 
define r 2 3 to be the distance where A23 = 1. This gives 



/ \ Lt,oi 
°( r ) = — 



1.14 x 10- n a(r) 8.8 x 10" 

( r / r23 )-2,/3 



r 



(ii) 



The column density dependent critical distance where 
radiation force dominates the cloud motion is, 



> 



7.7 x 10- 



a(r)T 



3/2s 



(12) 



For example, the case of s = 1 and a(r) = 0.5 gives a 
critical radius of r = 0.06r~ 15 7' 2 3 for radially moving 
clouds. This dependence of r on T is the main motiva- 
tion to suggest a new method for evaluating Mbh and 
tblr (§4). As explained, the critical radius should not 
be confused with the point of escape from the system. 
Non-radial velocity components (vg), that reflect the en- 
ergy and angular momentum of the system will act to 
reduce this radius (see examples below). 

The motion of BLR clouds with the above properties 
involves an acceleration term of the form, 



a(r) = 



cia{r) _ £2 

^•2 — 2s/ 3 y>2 ' 



where 



and 



ci = 1.14 x 10- n L fcoi r„ 2s/3 



c 2 = 8.8 x 10- 13 L bo; /r 
The radial potential is, 



$(r) 



a(r)dr 



(13) 

(14) 
(15) 

(16) 



where r* is the radius where $(r) = 0. Below we use this 
potential to calculate cloud orbits and line profiles. The 
energy and angular momentum terms that result from 
the above integration, are included in the calculation by 
fixing the initial conditions, r and the two velocity com- 
ponents at this location. 

3. LINE PROFILE CALCULATIONS 
3.1. Method 

We carried out a series of calculations under a variety 
of conditions considered to be typical of different BLRs. 
Every model is calculated for assumed Mbh and T. This 
specify Lj, i and thus the potential $(r). The additional 
model parameters are: 

1. The radial parameter s. 

2. The cloud column density normalization factor r 2 3. 

3. The initial radius ro and the initial velocity vq = 
vg(ro). We assume that the orbits of clouds with 
very large N co i are ellipses of given eccentricities. 
r is chosen to be the apogee of the orbit and 
vq (given below in units of the Keplerian velocity, 
v Kepler) is determined from these conditions. This 
a simple way to specify the angular momentum. In 
the examples below, we focus on those cases where 
the resulting FWHMs are consistent with the ob- 
servations of the broad H/3 and C iv A 1549 emission 
lines but give details for several others. 

4. The initial ionization parameter, U{tq). We note 
that the exact value of the gas density, rijj(r), is less 
important. In the following we assume n#(ro) = 
10 10 cm -2 for all cases. 

5. The three-dimensional distribution of orbits. This 
is done in two steps. First we calculate the motion 
of numerous identical clouds in a plane and then 
distribute many such planes in a spherical geom- 
etry specified by the inclinations of the planes to 
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Fig. 1. — Planar orbits of three clouds with V = 0.1 and different 
column densities. The large column density cloud (thick line, T2'i = 
lOOOro) moves in a closed elliptical orbit. A smaller column density 
cloud (thin line, r^'i = 10ro) moves in a closed rotating orbit and a 
marginal column density cloud (dashed line, T2Z = 0.82ro) escapes 
the system. 

the line of sight. The profiles given below are only 
those for a line of sight which is perpendicular to 
the central plane of motion (if such a plane exists). 
All calculations assume a large enough number of 
clouds such that the predicted profiles are smooth 
(see Bottorff and Ferland 2000 for discussion and 
earlier references on this issue). 

Physical properties that are not included in the present 
calculations are non-isotropic central radiation field, non- 
isotropic line emission, the photoionization of gas with a 
range of density and mctallicity, different inclinations of 
the line of sight to the central plane of motion, large 
cloud covering factors in a specific direction, and central 
obscuration, e.g. by the accretion disk. Several of those 
are likely to be important in real BLRs but are beyond 
the scope of the present work. 

Fig. [T] illustrates the orbits of three s — 1.2, tq = 10 17 
cm, r = 0.1 and vq = Q.bvKepier clouds moving under the 
influence of a 10 8 M Q BH. The first is an ellipse typical of 
a cloud which is not affected by radiation pressure (e.g. 
^23 = 1000ro). This is shown by a thick solid line. For- 
mally speaking, such clouds are Compton thick but this 
is of no practical implications since the only intention is 
to show a simple, gravity dominated orbit. The second is 
a case where = lOrrj. Here, radiation pressure force is 
significant and acts to constantly changing the direction 
of motion of the cloud. This results in a rotating planar 
orbit. The third orbit (dashed line) follows the trajectory 
of a smaller column density cloud (r 2 3 — 0.82ro) that es- 
capes the system. Increasing T will result in similar type 
orbits for the rotating orbit second cloud except that the 
angle between two successive revolutions will increase. 

We calculated various line profiles for the case of 
-^bh=10 8 M q , r 2 3 = 10r , vo — 0.5 and T in the range 
of 0.05 (negligible radiation pressure force) to 0.735 (just 
below escape). The bolometric luminosity in each of 
those is obtained from the combination of Mbh and T. 
We assume a(r) = 0.5 at all radii and e(r) which takes 
into account only geometrical factors (i.e. constant j(r), 
see eqn. \§§ and isotropic line emission. In terms of total 
line emission, this is a reasonable approximation for lines 
like H/3 that reprocess roughly a constant fraction of the 



Fig. 2. — Line profiles for spherical s=1.2 atmospheres around 
a 1O 8 M0 BH and a range of T as marked. All clouds start at 
ro = 10 17 cm with v r = and vg = 0.5 vjie-pier- The column 
densities are changing as {r/r23)~ 2s ^ 3 with T23 = 10 18 cm (N co i £3 
6.3-10 23 cm -2 at r ). The FWHM of the profile decreases with the 
increasing V due to the increasing importance of radiation pressure 
force. The profile parameters are listed in Table 1. 

ionizing continuum radiation. Obviously, a large opti- 
cal depth in H/3 will result in line emission anisotropy 
which is not considered here. It is not appropriate for 
lines like C iv A 1549 whose intensity is more sensitive to 
the level of ionization and the gas temperature. At this 
stage we specifically avoid the use of a varying a(r) since 
the effect on the orbit can be significant even for small 
changes in this parameter (see below). The resulting 
profiles, assuming a complete spherical atmosphere (the 
entire ±7r/2 radians range relative to the central plane), 
are shown in Fig. [2j As expected, the profile becomes 
narrower with the increasing T reflecting the fact that, 
as the luminosity increases, the cloud spend less and less 
time at small radii. 

The top part of Table 1 provides additional informa- 
tion about the calculations. For each profile we give 
the FWHM in units of VKepier(fo), the mean emissiv- 
ity weighted radius, < r > /ro, and the mass correction 
factor / (eqn. [2]). The calculation of < r > is obtained 
by weighting the emissivity of the cloud and the time it 
spends at each radius. This is roughly equivalent to the 
observed RM radius. The mass correction factor is ob- 
tained by requiring fFWHM 2 < r > /G = M B h- We 
also show (in parenthesis) the values of FWHM and / 
obtained for the case of a thick central disk which rep- 
resents only a part of a spherical distribution. Here the 
cloud distribution correspond to a width, relative to the 
central plane, of ±7r/4 radians. The reduction in FWHM 
relative to the complete sphere is about a factor of 0.6 
and there is a corresponding increase in /. Computed 
line profiles that are typical of this and similar geome- 
tries are shown in Fig. 01 

To explore models with different initial conditions, we 
computed two cases of planar orbits with the same or- 
bital energy and different angular momentum. One such 
example is shown in Fig. [3] The less eccentric case in 
the diagram corresponds to the orbit labeled with 2) in 
Fig. 1 for which vq = 0.5. The more eccentric one as- 
sumes vq = 0.25 but with a non-zero radial velocity of 
v r = 0.433. This results in a much narrower profile. In 
the middle part of Table 1, we report other cases where 
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TABLE 1 

Line widths, mass conversion factor /, and 
emissivity- weighted radii for various models 8, 



r 


FWHM/» Kefler (ro) 


< r > /ro 


/ 


s = 1.2 


T23 = 10r ( ) 


vq = 0.5 




0.05 


1.58 (0.93) 


0.54 


0.75 (2.18) 


0.1 


1.55 (0.92) 


0.54 


0.77 (2.21) 


0.3 


1.45 (0.87) 


0.56 


0.85 (2.37) 


0.5 


1.34 (0.81) 


0.59 


0.94 (2.56) 


0.7 


1.15 (0.72) 


0.68 


1.11 (2.78) 


0.735 


1.06 (0.68) 


0.78 


1.13 (2.76) 


s = 1.2 


r 2 3 = 10r 


Vq = 0.25 




0.05 


1.04 


0.45 


2.05 


0.1 


1.02 


0.45 


2.10 


0.3 


0.95 


0.47 


2.39 


0.5 


0.87 


0.49 


2.74 


0.7 


0.76 


0.52 


3.31 


0.91 


0.59 


0.67 


4.32 


s = 1.2 


r 2 3 = ro 


vo = 0.5 




0.01 


1.57 


0.54 


0.76 


0.03 


1.51 


0.55 


0.80 


0.1 


1.23 


0.64 


1.03 


0.116 


1.055 


0.79 


1.13 



a Assuming the line emissivity is strictly proportional to the cloud 
cross section and a(r) = 0.5. In all cases vo = vg(ro). Numbers 
for / assume spherical BLRs (numbers in brackets assume a ±7r/4 
radians thick disk). 



vq = 0.25 and v r = 0. Such orbits are again very ec- 
centric and the profiles are, indeed, much narrower. The 
corresponding values of / are now larger by a factor of 
2-3 than those observed. Additional models (not shown 
here) with larger initial angular momentum, give larger 
FWHM and smaller /. Obviously, some combination of 
all those is required to explain real observations. The 
bottom part of Table 1 shows the results of a set of line 
profile calculations carried out for smaller column den- 
sity clouds. We chose r-zi = r o which corresponds to a 
factor of 6.3 decrease in iV co i relative to the case shown 
at the top of the table. The scaling of FWHM between 
the two cases is simply by the corresponding factor in T 
(i.e. the same FWHM for V smaller by a factor of 6.3). 
This illustrates the fact that in an atmosphere with a 
large range of column densities, there are always clouds 
that are close to being ejected from the system at large 
distances. 

The changes in < r > for a given r 2 3 shown in Table 1 
are due to the fact that as T increases, and radiation 
pressure is more important, the clouds spend more and 
more time away from the BH. This is noticeable for the 
case of r23 = 10ro when T approaches 0.73 and for the 
case of r23 = ro when T approaches 0.1. 
As noted in §1, RM campaigns show that rB£^(H/3)oc 



L 



0.6±0.1 
bol 



It is interesting to note that this behavior is 
not very different from what is calculated here for the 
changes in < r > if we compare values over the range 
where T approaches its limiting value. However, it is not 
the case when T changes by similar factors close to the 
lower range shown in the table, where radiation pressure 
is negligible. 

The values of / computed here should be compared 
with those determined observationally for selected AGN 
samples with measured cr*, in particular the Onken et 
al. (2004) and the Woo et al. (2010) AGN samples. 
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Fig. 3. — Top: Planar orbits of two clouds with the same orbital 
energy and different angular momentum. The model parameters 
are T = 0.1, s = 1.2 Mbh= 10 8 Mq and other parameters as in 
case 2) of Fig. 1. The less eccentric case ((1), thin line) corresponds 
to maximum angular momentum at r = ro with do = 0.5 and no 
initial radial velocity (tv,o = 0). The more eccentric case ((2), 
thick line) assumes vq = 0.25 and t> r ,o = 0A33vk epler ■ Bottom: 
Line profiles for the two cases (same notation as in Fig. 2). The 
narrower profile corresponds to orbit 2). 



The simulations illustrate how this factor depends on 
the BLR geometry, the distribution of T among objects 
in the sample and the distribution of A col in individual 
BLRs. 

An important point of the new calculation is the rela- 
tively little change in the value of FWHM listed in Ta- 
ble 1, only a factor of <~ 1.5 over most of the range of 
T except very close to the limiting value. The changes 
in / are also small, only a factor of ~ 1.3 over the same 
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Fig. 4. — Same initial conditions as in Fig. [2] for T = 0.5. The 
various profiles represent motion in different spherically shaped at- 
mospheres. The narrowest profile (dashed line) represents a sphere 
where clouds occupy only the section between -0.3 and +0.3 radi- 
ans relative to the mid-plane (which is perpendicular to the line of 
sight). The other cases are for wider coverage with clouds between 
-0.9 and +0.9 rad (dotted line) and -L5 to +L5 rad (solid line). 
The double peak profile illustrates the case of two polar caps where 
the clouds occupy a sphere whose mid-section, between -L2 and 
+1.2 rad, has been removed. 

range in T. This seems to be in contradiction to the 
naive expectation that, for cases of increasing L, the 
term < r > FWHM 2 /G will deviate more and more 
from A/bh (e.g. eqn. [3]). There are two reasons for this 
behavior. First, for realistic cases where iV co i depends 
on the cloud location, the mean emissivity distance and 
the velocity depend on L/M rather than on L. This 
suggests that very low and very high luminosity AGNs 
with similar T will react to radiation pressure force in a 
similar way. Second, for a planar motion, the changes in 
the radial potential <fr(r) do not affect the cloud velocity 
in a linear way. In fact, the mean orbital changes in vg 
are small enough such that the overall FWHM is very far 
from zero even for marginally escaping clouds. Moreover, 
the mean cloud location, < r >, is increasing in reaction 
to the increasing radiation pressure term. The end re- 
sults is that the product f < r > FWHM 2 /G, with 
a constant value of /, is always a reasonable approxi- 
mation for Mbh with little dependence on the relative 
importance of gravity and radiation pressure force. We 
return to this issue in §4 where we suggest a new way 
to evaluate Mbh taking into account radiation pressure 
acceleration. 

Finally, we note that while radiation pressure is neg- 
ligible for very small values of T, the s-dependence of 
the cloud properties is still very important. For exam- 
ple, an s = atmosphere gives constant column density 
clouds (similar to what was assumed in M08) yet, the 
mean emissivity radius, the FWHM of the emission lines 
and the mass correction factor / in this case are always 
different from those of the s = 1.2 case, regardless of 
the column density. The reason is the dependence of the 
cloud cross section on s. For example, in the case of 
T = 0.01 (first entry in the bottom part of Table 1), the 
s = case gives < r > /ro = 0.39 (compared with 0.53 
for s = 1.2) and FWHM/vKepier = 2.45 (compared with 
1.51). The resulting / is therefore much smaller (0.42 
compared with 0.76). Thus, the radial dependence of 
the cloud properties are important for all T. 



3.2. Applications to spectroscopic observations of 
AGNs 

The examples discussed above were normalized to give 
a typical r BLR (R(3) for AGNs with M B h=1O 8 M and 
T = 0.1. However, the computed line profiles cannot be 
directly compared with the observations of such sources 
for several reasons. First, we only consider a situation 
involving one type of clouds and neglect the possibil- 
ity of different populations under different physical con- 
ditions in the same source. This applies to the distri- 
butions of both A co i and U(r). For example, eqn. Q] 
and the constants given in §1 suggest that, in gen- 
eral, rBLR(Hp)/rBLR{Civ\15AQ) ~ 3. The question is 
whether cloud distributions like those considered in Ta- 
ble 1 can reproduce this ratio. Second, we did not take 
into account changes in a(r), the fraction of Lf, i which 
is absorbed by clouds at various distances. This can be 
an important factor close to the BH where clouds be- 
come partly transparent. In this case, much of the Ly- 
man continuum radiation is not absorbed and radiation 
pressure force is reduced. It can also affect medium to 
large column density clouds at large distances where a(r) 
approaches unity. For example, assuming a{r) — 0.75 in- 
stead of a (r) — 0.5 in the calculations of Table 1 results 
in a limiting value of Y which is about 0.4 compared with 
T = 0.735 listed in the table. 

To illustrate these effects, and to provide more realistic 
line profiles, we computed two grids of photoionization 
models for a range of column density and ionization pa- 
rameter using the code ION (Netzer 2006). The first grid 
supplies calculated line intensities for H/3 and CivA1549 
over a large range in U(r). Given r from the cloud mo- 
tion simulation, we use the grid to compute j(r) and 
thus a more realistic e(r). The second grid supplies the 
absorbed fraction, a(r), as a function of U(r) and iV co i. 
Fig. [5] shows part of the a(r) grid to illustrate the ex- 
pected range in this parameter. We have not included 
the changes in gas density since they do not play a ma- 
jor role over the range of conditions considered here. We 
have also not considered anisotropy of the line emission 
which is bound to have an effect on the FWHM of some 
lines. Such modifications will be included in a forthcom- 
ing paper that is intended to present a comparison with 
observed line profiles. 

We tested a large number of single-zone models 
using the above grids of U(r) and a(r). The models 
cover a large range in angular momentum and BLR 
geometries. We have specifically investigated three 
cases of different eccentricity, defined by three values 
of vo(ro), 0.25, 0.5 and 0.75. These were calculated 
with different Y and r^-i- In general, it is easy to 
reproduce the observed I(ClV A1549)/I(H/3) but diffi- 
cult to account, at the same time, for the emissivity 
weighted radii of the two lines (eqn. [1} and the line 
width ratio. For example, the best case of the three, 
with v (r ) = 0.25, gives I(C IV A1549)/I(H/3)=4.4, 
< r >(CivA1549)/< r >(H/3)=0.67 and 
FWHM(CivA1549)/FWHM(H/3)=1.43. The con- 
clusion is that, within the range of parameters assumed 
here, there is no obvious way to explain all those 
properties when keeping with the idea of a single column 
density distribution (i.e. a single ^3). 

We also tested a case of Mbh= 1O 8 M , Y = 0.1 and 
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Fig. 5. — Part of the a(r) grid (fraction of the total contin- 
uum flux absorbed by the clouds) used in the present calculations. 
Numbers along the contour lines are log a. 



TABLE 2 

Properties of the two-zone model with uo(ro) = 0.5. 



Zone FWHM(H/3)FWHM(Cto) 
(kms -1 ) (kms -1 ) 


r(U0) 


r((Jiv ) 


l(Civ) 


r(RM,H0) 


r(RM,Civ) 


J(H/3) 


Inner 3160 


3450 


0.32 


0.88 


9.55 


Outer 1390 


2580 


1.64 


3.2 


1.45 


Combined 2060 


3390 


0.98 


1.1 


5.5 



two distinct cloud populations in inner and outer zones 
with some overlap between the two. In this case, the 
initial conditions for the two populations are decoupled 
from each other but the changes in density, column den- 
sity and ionization parameter follow the same pattern 
with the same s = 1.2 density law. The FWHM of 
both emission lines were calculated under the assump- 
tion of a thick central spherical sector with clouds oc- 
cupying a region of ±7r/4 radians relative to the central 
plane. The inner zone clouds have r$ = 5 x 10 16 cm and 
U(ro) = 10 _1 and the outer-zone clouds tq = 3 x 10 17 cm 
and U(ro) — 10~ 25 . The starting velocity in both zones 
is vq = 0.5 at the appropriate r$. In both zones r23 = 3ro. 
We followed the cloud motion and calculated, in each 
zone, the line intensity ratio, I(C IV A1549)/I(H/3), the 
line FWHMs, and the emissivity weighted radii. These 
numbers are listed in Table 2 where we also show the 
properties of the combined spectrum which is calculated 
under the assumption of equal contributions to H/3 from 
both zones. The emissivity weighted radii for the two 
zones are given in units of the RM-radii of the two lines 
(cqn. [1} . This very simple two- zone model gives results 
that are in good agreement with the observations of many 
low-to- intermediate luminosity AGNs. Fig. [B] is a graphi- 
cal summary of these results. The left and central panels 
show H/3 and CivA1549 profiles for the inner and outer 
zones, again assuming isotropic line emission, and the 
right panel shows the combined two-zone profiles. 

In conclusion, the simple single zone models explored 
here cannot reproduce all the observed properties: line 



intensity ratio, mean emissivity radii and FWHM ratio. 
The main reason is that the starting conditions fix the 
cloud orbit, and hence the line emissivity and FWHM. 
Simple two-zone models like the ones presented here 
can account for most observed properties of the H/3 and 
Civ A 1549 lines. In particular, they can account for the 
mean line ratio, the mean emissivity weighted radii and 
the mean relative FWHM of the H/3 and C iv A1549 lines 
measured in various RM samples. Obviously, such sim- 
ple models do not intend to explain all the observed line 
profile properties that can differ from one object to the 
next and contain additional components (see some obvi- 
ous examples for complex C IV A 1549 profiles in Richards 
et al. 2002 and Sulentic et al. 2007). Fitting those is 
deferred to a forthcoming paper. 

4. DISCUSSION 
4.1. General considerations 

The above calculations allow us to investigate the in- 
tensity, the width and the shape of the broad emission 
lines and to evaluate various methods used to estimate 
Mbh- We defer the discussion of specific observed line 
profiles to a future paper. 

Assume a system of clouds with a given total amount of 
gas and a large range of column densities. Such a system 
will eventually break into three: virialized clouds, non- 
virialized bound clouds and escaping clouds. The third 
group will not contribute significantly to the observed 
line emission for more than several dynamical times. The 
relative contribution of the first and second groups to the 
line emission depend on the cloud mass distribution. A 
sudden increase in Lb i will increase the importance of 
radiation pressure and will remove more gas from the 
system. A decrease in L^oi will drive the system closer 
to virial equilibrium. A new gas supply, e.g. from a disk- 
wind, will produce bound as well as unbound clouds. 
All aspects of this general scenario must be considered 
when evaluating the observed line profiles and the various 
methods developed to use them in estimating Mbh- 

A major objective of the present paper is to evaluate 
the accuracy and the normalization of various Mbh esti- 
mators in type-I AGNs. The results presented in Tables 
1 & 2 suggest the following: 



Netzer, Marziani 



1. Every AGN is likely to contain a large number of 
clouds with a large range in N co \. This can be the 
result of a broad cloud mass distribution and/or 
due to cloud motion in a radial-pressure dependent 
environment with a positive value of s. A given T 
results in a lower limit on iV C oi at a given location 
for a given orbit eccentricity. Under such condi- 
tions, there are always some clouds, e.g. those that 
are very close to the BH, for which radiation pres- 
sure is negligible. For others, radiation pressure 
can be very important. 

2. For a small enough iV co i, the effective tblr depends 
on both T and iV co i. Under these conditions, the 
BH mass itself is an important factor in determin- 
ing tblr- To illustrate this, consider two AGNs 
with identical SED, Lboi, BLR geometry, Ar co i dis- 
tribution and inclination to the line of sight. The 
effective tblr. in the two is the same provided they 
harbor identical BHs. Different tblr will be mea- 
sured if the two BHs have different masses despite 
of the fact that Li, i is the same in both. This is the 
result of the larger T in the smaller BH AGN. The 
effect may not be recognized in a large sample of 
sources and can, in fact, be attributed to a large in- 
trinsic scatter in the Lb a i — tblr relationship. Any 
derived L^oI — tblr relationship will depend on the 
properties of the sources in the chosen RM sample, 
in particular on the distribution of T. 

3. Assuming a range in iV co i in every AGN, the M08 
suggestion to include a luminosity dependent term 
in the calculation of Mbh (eqn. [3]) is not in accord 
with our calculation that indicate that tblr and 
FWHM depend on L/M and not on L. 

The multi-year RM campaign of NGC 5548 is the best 
example to test some of these ideas in a specific source. 
The campaign has been described and analyzed in nu- 
merous papers and the ones most relevant to the present 
study are Peterson et al. (1999) and Gilbert and Peter- 
son (2003). 

Fig. [7] shows the variations in L5100 and time lag (in 
this case the centroid of the CCF) in NGC 5548. Each 
point represents a full observing season which is typi- 
cally ~ 300 days long. The data are taken from the 
recent compilation of Bentz et al. (2009) which provides 
the best galaxy subtracted flux at 5100A. The uncer- 
tainty on L5100 is basically the range of this quantity 
over the observing season. This is of the same order as 
the variation from one season to the next. As clearly seen 
from the diagram, rsLR^P) lags the continuum in such 
a way that more luminous phases are associated with 
longer lags. This has been noted in earlier publications, 
e.g. Gilbert and Peterson (2003). Fig. [8]shows t(lag) vs. 
£5100 for the same data set. While the uncertainties are 
large, some correlation, with a slope of 0.5-1, is evident. 
An earlier version of the diagram, with fewer points, is 
shown in Peterson et al. (1999). 

For NGC 5548, Mbh— 10 8 M© and r BLR (UP)~ 20 
l.d.. Thus, the dynamical time is of order 6 years and 
the time it takes to change tblr by 50% (e.g. Table 1) 
is approximately 3 years. This seems to be compatible 
with the changes in L5100 and t(lag) in fig. [3 thus some 
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Fig. 7. — Changes in continuum luminosity (-L5100) an d time lag 
for NGC 5548 (data from Bentz et al. 2009). Error bars on L5100 
were omitted for clarity. 
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Fig. 8. — The correlation of L5100 vs. t(lag) for NGC 5548. Data 
is in fig. The dashed line has a slope of 0.5. 



adjustment of rsLR^P) due to the effects discussed in 
this work are possible. The measured £5100; with a bolo- 
metric correction factor of about 10, indicates a mean 
r of about 0.02. The bottom part of Table 1 provides 
approximate parameters for such a case. Any success- 
ful model of NGC 5548 should account for the behavior 
shown in Fig. [3 as well as for the observed FWHMs and 
luminosities of both H/3 and CivA1549. While the full 
investigation is deferred to a future paper, we consider 
here the predicted lags for T =0.005, 0.01, 0.02 assum- 
ing a BH mass of 10 s M Q and two families of clouds: one 
with r23 = 10ro and vq = 0.5, and one with r23 = 0.093ro 
and vq = 0.25. The second assumed family of clouds re- 
sults in pseudo-orbits of higher eccentricity that, as ex- 
plained earlier, are more strongly affected by radiation 
pressure. In both cases, the predicted lags for T — 0.02 
are consistent with the observed values (logi(lag) « 1.34 
at logisioo ~ 43.48). However, the calculated slope of 
logi(lag) vs. L5100 is flatter than observed. 

As argued earlier, a single family of BLR clouds can- 
not provide a full explanation to the observed spectrum 
of many AGNs. This must applied to NGC 5548 (to 
appreciate the complexity of this case see the various 
components considered by Kaspi and Netzer (1999) to 
explain only the variable line intensities). The simple 
examples considered here suggests that dynamical seal- 
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ing of the BLR in NGC 5548, due to radiation pressure 
force, is an additional, physically-motivated mechanism 
that must be added to any cloud model when attempting 
to explain the observed variations in i(lag). 

4.2. Evaluation of present Mbh estimators 

Current BH mass estimates utilize RM-based measure- 
ments of rBLR, measured FWHMs (or an equivalent ve- 
locity estimator) of certain broad emission lines, and 
cqn. [2j The normalization constant / is obtained by 
a comparing Mbh obtained in this way with the mass 
obtained from the M — c* method. Having examined a 
large range of cloud orbits and line profiles under various 
conditions, and the corresponding values of the effective 
tblr, we are now in a position to evaluate the merits of 
this method. 

We consider three general possibilities. The first is the 
case where all AGNs contain BLR clouds with a wide 
column density distribution. A randomly chosen object 
will have in its BLR some clouds that arc affected by 
radiation pressure force and others that are not. This 
is the case for any F. The cloud dynamics and the ob- 
served line profiles reflect the (unknown) column density 
distribution. Our calculations suggest that an RM sam- 
ple drawn randomly from such an AGN population can 
be safely used to determine the best value of / by com- 
paring the derived Mbh with the M— <j* method. This is 
justified by the fact that Mbh °c< r > FWHM 2 even if 
radiation pressure force is important (see Table 1). The 
observed FWHMs are, indeed, smaller than the ones that 
would have been observed if all clouds had extremely 
large column densities. This, however, has no practical 
implication since the column densities are not known and 
/ is simply a normalization factor that serves to bring 
two completely different methods of estimating A/bh into 
agreement. Mass estimates obtained in this ways are reli- 
able provided the properties of the RM sample represent 
well the population properties. 

The second case reflects a situation where the cloud 
column density distribution is, again, very broad but part 
of the population is under-represented in the RM sample. 
For example, the RM sample may contain mostly sources 
with r ~ 0.1 while the overall distribution of L is much 
wider. In this case, the normalization factor / will reflect 
only the properties of the measured sources and its use 
will provide poor mass estimates for cases with much 
larger or much smaller accretion rates. This may well 
be the case in the RM sample which is most commonly 
used (Kaspi et al. 2000; Bentz et al, 2009) that contains 
only very few AGNs with T > 0.3. The numbers in 
Table 1 enable us to evaluate the resulting deviations 
in the estimated Mbh- For example, if we use the first 
part of the table and assume a source with a certain L\, Q i 
and r = 0.1, we find that the mass of a similar Lf, i 
source with L = 0.7 will be under-estimated by a factor 
of 1.11/0.75. 

Regarding the second case, it is important to note that 
under-estimates and over-estimates of Mbh are equally 
likely. Consider again an RM sample where, for most 
sources, T = 0.1. This results in a certain value of / 
which takes into account the effect of radiation pressure 
force in some of these sources (see bottom part of Ta- 
ble 1). Assume a second, randomly selected AGN sam- 
ple with a similar BH mass distribution but a typical L 



which is much smaller than 0.1. Most measured FWHMs 
in this sample are broader than those in the RM sample 
because radiation pressure force is not as effective in re- 
ducing the cloud velocity. Using the value of / derived 
for the RM sample will result in over estimating Mbh 
in the second sample. The lower part of Table 1 gives 
some idea about the magnitude of this effect, e.g. an 
over estimate by a factor of 1.01/0.76. 

The third case is similar to the first one except that 
large luminosity variations, on time scales that are not 
too different from the BLR dynamical time, are occur- 
ring in most sources, including those selected for RM 
monitoring. Table 1 shows that, like the first case, the 
deduced / represents well the population because < r > 
follows the variations in L{, i. The mean Mbh in such a 
sample is recovered albeit with a larger uncertainty. 

4.3. Alternative Mbh estimators 

Given the above considerations, we now investigate an 
alternative method to calculate Mbh- The method takes 
into account the effect of radiation pressure force on the 
cloud motion and the results will be compared to those 
obtained by the old method (eqn. [5]) and by the M08 
method. 

Our new calculations indicate that the emissivity 
weighted tblr depends both on the (large) range in L 
across the entire AGN population, as well as on short 
time scale changes in rsLR in individual sources. The 
first of those depends roughly on L 1 / 2 and is a manifesta- 
tion of the observational fact that the ionization parame- 
ter, U(r), and the spectral energy distribution (SED), are 
not changing much with source luminosity. The second 
reflects changes in the BLR structure due to the reac- 
tion of various column density clouds to the (changing) 
radiation pressure force. This depends on both Lt> ; and 
Mbh- This is seen for example in eqn. 1121 for the critical 
radius where clouds can escape the system and also in 
the calculations of Table 1. It is therefore reasonable to 
assume that rsLR is given by an expression of the form, 

rsLR^a^ + a 2 (L/M) s , (17) 

where a\ and ai are constants and L is a measure of 
the source luminosity, e.g. L 5 i o if Tblr = Tblr (H/3) . 
Obviously, the above approximation is not unique and 
one can assume other dependences that are consistent 
with the line profile calculations, e.g. a dependence of 
FWHM on L/M. 

The idea of introducing a second, luminosity depen- 
dent term into the calculation of Mbh is not new. In 
particular, M08 suggested an expression for Mbh which 
depends on both L 1 / 2 and L/N co i (eqn. [3]). Assuming all 
AGNs obey the same relationship, and JV C oi is the same 
in all, the M08 expression leads to extremely large values 
of Mbh for the most luminous AGNs. The reason is the 
linear dependence of Mbh on L at very high luminosities 
combined with the calibration of the relationship at small 
L, typical of the M — ex* sample of Onkcn et al. (2004). 
The additional consequence of this approach is an upper 
limit of r ~ 0.1 in many high luminosity, large BH mass 
sources. In their later work, M09 considered the possi- 
bility that JV co i can differ from one source to another but 
is still constant for all clouds in a given BLR. This would 
result in smaller Mbh and larger T in some high lumi- 
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nosity sources since in some BLRs, N co \ can exceed 10 23 
cm -2 by a large factor thus reducing the importance of 
radiation pressure force. 

The limitation of the M08 mass estimate is the detach- 
ment of L from A/bh- As shown here, this is not the case 
in more realistic BLRs, especially those where the masses 
of the clouds are conserved. In such cases, the location of 
the outer clouds that still contribute to the line profiles 
depends on L/M and the 3D-velocities of the marginally 
bound clouds are such that the product tblrFWHM 2 
is not very different from what is found in pure grav- 
ity dominated systems. Moreover, for pressure confined 
clouds, the dependence on N co \ is likely to be different 
in different parts of the BLR. Thus, we are looking for 
an expression that will reflect, properly, all these effects 
and will allow for the possibility of a range of column 
densities in every source. We also want to avoid biasing 
in the derivation of Mbh in the limits of very large or 
very small L and to retain the experimental results that 
tblr oc L 1 with 7 = 0.6 ± 0.1. 

All the above can be achieved by assuming that 
tblr is given by eqn. [T7] and requiring that Mbh oc 
rBLRFWHM 2 . For the sake of simplicity, we assume 
7 = 0.5 and (5 = 1 and substitute eqn. [17] into the mass 
expression. This leads to a simple quadratic equation in 
Mbh with the following solution, 



M 



BH 



= -a 1 L 1 l 2 FWHM' 1 
2 



n 2 \ 

1 + V(!+ a 2 FWHM 2) 

(18) 

where a\ and a 2 are the same ones used in eqn. [171 except 
for a common multiplicative constant which depends on 
the units of rBLR, L5100 and Mbh- For example, us- 
ing the measure parameters for the H/3 line, L =£5100, 
FWHM=FWHM(H/5), then the constant multiplying ai 
and a 2 in eqn. [TTJis 10 16123 when M B h is measured in 
Mq, L5100 in units of 10 44 ergs s~ x and rsLR in cm. 

We used eqn. [18] and the Woo et al. (2010) sample to 
find cti and a 2 for 29 AGNs with measured <r*. The list is 
an extension of the one used by Onken et al. (2004) that 
contains only 16 sources. We have supplemented the data 
in Woo et al. by data from Bcntz et al. (2009) on r B LR 
and L5100 where this information was missing. First, we 
performed a x 2 analysis on Mbh(RM) vs. M — a* using 
the parameters recommended by Giiltekin et al. (2009). 
This gave / = 1.0 which is consistent with the values 
found by Onken et al. (2004) and Woo et al. (2010)0 

Next we carried out a x 2 minimization to solve for 
ai and a 2 in eqn. 1181 Since the minimization involves 
the error estimate on L5100, and since this error is not 
very well defined given the combination of observational 
uncertainly and the intrinsic scatter in L5100 over sev- 
eral long RM campaigns, we decided to adopt a uniform 
value of AL5ioo/-t'5ioo = 0.3. We also assume a min- 
imum of 0.1 to A(cr*)/er* and a minimum of 0.05 on 
A(FW H M) / FW H M . Our results depend slightly on 
these assumptions. 

3 Onken et al. (2004) and Woo et al. (2010) carried the anal- 
ysis using the H/3 line dispersion rather than FWHM(H/3). For 
the sample in question, this line- width measure is smaller than the 
FWHM(H/3) by a factor of approximately 1.9 leading to a corre- 
sponding increase in the mean / by a factor of about 1.9 2 . All 
these numbers are sensitive to the error estimate in cr* and in the 
virial product. 



The best values obtained in this procedure are a\ 
■1.1. a 2 = 7.1 x 10 7 and x 2 l v — 1-73. Extensive tests show 
that the x 2 changes very little if a\ or a 2 are changing by 
up to 10%. This is the result of some degeneracy between 
cii and a 2 (see eqn. II 8|) . The average deviation between 
the new mass estimates and those obtained by the A/— cr* 
method is 0.31 dex. There is a weak dependence of the 
deviation on the line width (larger deviation for larger 
FWHM(H/3)) which is marginal given the small number 
of sources in the sample. The corresponding number for 
the deviation of masses obtained directly from the RM 
measurements and the above value of / is 0.36 dex. Thus 
the new method is, indeed, superior in this respect. Ob- 
viously it is not surprising to find such an improvement 
when adding a new free parameter to the model. 

To compare the various mass estimates more thor- 
oughly, we calculated Mbh in three different ways: the 
old method (eqn. [2]) with / = 1.0, the M08 method 
(eqn. \5§ with /i = 0.56 and f g = 10 77 (as in M08), 
and the new method feqn.fTSj) with the above a\ and a 2 . 
For the M08 method, we followed the M09 recommenda- 
tion and assumed a log normal distribution of N co \ with 
a mean of 10 23 cm -2 and a large standard deviation of 
0.5 dex. We also calculated rsLR in the old (eqn. [T]) and 
new (eqn. I17[) ways. 

Fig.|9]compares two mass ratios, Mbh (new) /Mbh (old) 
(red points) and Mbh (new)/MsH (M08) (black points), 
in a large simulated AGN sample. The sample cov- 
ers, uniformly, the luminosity range Lsioo^lO 43 — 10 47 
ergs s _1 and the simulations assume a Gaussian, lumi- 
nosity independent distribution of FWHM(H/3) with a 
mean of 4,500 km s _1 and a variance of 1,500 km s -1 . 
The diagram shows that the new and old estimates are 
similar at all Lf, i but Mbh(M08) deviates from both, by 
a large factor, at both low and high luminosities. More- 
over, the slight deviation between the new and old meth- 
ods at the very high luminosity end, by up to about 0.2 
dex in Mbh, is most likely due to the fact that the proce- 
dure used to obtain a\ and a 2 is based on a sample of 29 
mostly low-to-intcrmediate luminosity AGNs while the 
simulations reach a much larger value of L\, D i. A com- 
parison of the estimated rsLR (eq. Q] and QTJ leads to 
similar conclusions. 

We also made a similar test on the Netzer and Trakht- 
enbrot (2007) sample using all three methods. The lumi- 
nosity range in this case is smaller but the FWHM(H/3) 
distribution more typical of observed AGNs. The results 
(not shown here) are very similar to those of the simu- 
lated sample. 

In conclusion, the new method for estimating Mbh 
gives results that do not deviate much from the old 
method which is based on a single constant /. This is 
true at both high and low luminosities and over a large 
range in FWHM. Obviously, the range of parameters 
tested here (s, orbit eccentricity, several types of cloud 
distributions, etc.) is rather limited and more extensive 
modeling is required to confirm these results. However, 
it is our opinion that the main limitation of the Mbh 
determination methods remains observational and is re- 
lated to the fact that the present AGN M — a* sample 
is small (29 sources) and cannot possibly represent the 
entire range of properties, mostly T, observed in AGNs. 
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Fig. 9. — Comparison of the various methods for calculating 
A/bh- The ratio of the new-to-old (red) and new-to-M08 (black) 
methods are shown as a function of Lj, D ; for the simulated sample 
described in the text. Note the good agreement between the old 
and the new methods and the large deviation from the method 
described in M08 for very large and very small values of L bo i. 

5. CONCLUSIONS 

We have investigated the motion of BLR clouds with 
time-independent mass under a range of conditions de- 
fined by a radial-dependent confining pressure. These 
conditions enforce a range of A^ co i in every BLR, even 
if the intrinsic mass distribution of the cloud is nar- 
row. We calculated cloud orbits under a central potential 
that includes a radiation pressure term. The orbits were 
then combined to predict emission line profiles in several 
simple situations. We only considered uniformly emit- 
ted emission lines and the preliminary comparison with 
with actual observations used realistic emissivity and col- 
umn density distributions but was limited to the H/3 and 



CivA1549 lines and at most two different cloud distri- 
butions. We found significant changes in cloud locations 
and velocities for those cases where the column densities 
are small enough to allow a significant contribution due 
to radiation pressure. This can be important in both 
high and low T sources. However, while cloud orbits are 
strongly influenced by the radiation pressure force, there 
is a relatively small change in the mean tblrFWHM 2 
and hence no large underestimation or overestimation of 
Mbh- We illustrate this behavior in several cases but 
note that other cloud distributions, with different mass, 
location and velocity distributions, may lead to some- 
what different conclusions. We used the new results to 
suggest a novel method for calculating tblr and Mbh 
by applying two new constants that were calculated by 
a comparison of the H/3 and L5100 observations and the 
M - a* AGN sample of Woo et al. (2010). We applied 
the method to several large observed and simulated AGN 
samples and demonstrated good agreement between the 
new and the old, pure gravity based methods. The com- 
parison with the M08 methods shows large deviations in 
the estimates of A/bh • 
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